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1 Introduction 

The concept of Lie algebroid was introduced by Pradines in [231 [24] and has proved to be 
a useful tool in the formulation and analysis of many problems in differential geometry and 
applied mathematics |13[ [2]. In the context of geometric Mechanics, a program was proposed 
by A. Weinstein [25] in order to develop a theory of Lagrangian and Hamiltonian systems on 
Lie algebroids and their discrete analogs on Lie groupoids. In the last years, this program has 
been actively developed by many authors, and as a result, a powerful mathematical structure is 
emerging. The purpose of this paper is to review some of such recent developments. 

One of the main features of the Lie algebroid framework is its inclusive nature. In what 
respect to Mechanics, under the same formalism one can describe such disparate situations as 
Lagrangian systems with symmetry, systems evolving on Lie algebras and semidirect products, 
or systems with holonomic constraints (see |12[ [9] for recent reviews) obtaining in such cases 
Lagrange-Poincare equations, Poincare equations, Euler-Poincare equations or Euler-Lagrange 
equations for holonomically constrained problems (see [5j [6] ) . 

While the Lie algebroid approach to Mechanics builds on the geometrical structure of the 
prolongation of a Lie algebroid [15] , the origin of Lagrangian Mechanics is the calculus of varia- 
tions. It is therefore important to have a variational description of Lagrange's equations for 
a Lagrangian system defined on a more general Lie algebroid. We will show that Lagrange's 
equations for a Lagrangian system on a Lie algebroid are precisely the equations for the critical 
points of the action functional defined on the set of admissible curves on a Lie algebroid with 
fixed base endpoints, and we will also show how to find such equations by means of a Lagrange 
multiplier method [20] . 

One of the advantages of such a unifying formalism it that morphisms establish relations 
between these apparently different systems, leading to an adequate way to study reduction 

*This paper is a contribution to the Proceedings of the Workshop on Geometric Aspects of Integ- 
rable Systems (July 17-19, 2006, University of Coimbra, Portugal). The full collection is available at 
|http: / / www.emis.de/journals / SIGMA/Coimbra2006.html | 
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theory. In particular we will show how to reduce the variational principle and the symplectic 
equations in presence of a fiberwise surjective morphism of Lie algebroids. 

The extension of this ideas to the theory of optimal control systems was initiated in |17] . 
and will also be briefly reviewed. On any Lie algebroid a generalized version of Pontryagin 
maximum principle can be established in a global and coordinate free way which stresses its 
geometric properties and can be successfully reduced under morphisms. 

There are many other interesting aspects about the application of Lie algebroids to Mechanics 
which are not covered in this review. For other applications to control theory [7J, to discrete 
mechanics [TJ] and field theory [THJ, [TO] see the recent review [§]. For extensions to time- 
dependent mechanics see [2T] , 

The paper is organized as follows. In Section [2] we present some basic facts on Lie algebroids, 
including results from differential calculus, morphisms and prolongations of Lie algebroids. In 
Section [3j we give a brief review of the Hamiltonian and the Lagrangian formalism of Mechanics 
on Lie algebroids. In Section [J] we show that Lagrange's equation for a Lagrangian system on 
a Lie algebroid can be obtained by means of variational calculus by selecting an appropriate 
class of variations. Much inside is gained by studying the geometry of the infinite dimensional 
manifold of admissible curves, which is done in Section [5j In Section [6] we study the transfor- 
mation rules induced by morphism of Lie algebroids on the geometric objects of the theory, and 
how this is useful in the theory of reduction of Lagrangian systems. Finally in Section [7J we 
show how Pontryagin maximum principle can be extended for control systems defined on Lie 
algebroids and how to reduce optimal control problems. 

2 Preliminaries 

Lie algebroids 

A Lie algebroid structure on a vector bundle r : E — > M is given by a vector bundle map p: E — > 
TM over the identity in M, called the anchor, together with a Lie algebra structure on the 
C°°(M )-module of sections of E such that the compatibility condition [a, frj\ = (p(cr)f)rf+f[a, rj] 
is satisfied for every / € C°°(M) and every a,rj £ Sec(E). See [21 [13] for more information on 
Lie algebroids. 

In what concerns to Mechanics, it is convenient to think of a Lie algebroid as a generalization 
of the tangent bundle of M. One regards an element a of £ as a generalized velocity, and the 
actual velocity v is obtained when applying the anchor to a, i.e., v = p(a). A curve a : [to,ti] — > E 
is said to be admissible or an E-p&th if ^(t) = p(a(t)), where 7(t) = r(a(t)) is the base curve. 

A local coordinate system (x*) in the base manifold M and a local basis {e a } of sections 
of E, determine a local coordinate system (x l ,y a ) on E. The anchor and the bracket are locally 
determined by the local functions p % a and CSL, on M given by 



The functions p l a and Cg satisfy some relations due to the compatibility condition and the 
Jacobi identity which are called the structure equations: 



and [e a ,ep] = e r 



Pi 



Pa 




+ C^ 7 Ca/i + OyaCpn + C^C^ — 0. 
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Cartan calculus 

The Lie algebroid structure is equivalent to the existence of an exterior differential operator 
on E, d: Sec(A k E*) -» Sec(A fc+1 £*), defined as follows 

k 

du(a , . . . , a k ) = ^(-l)V(o-»)(w(cro, • • • , <H, ■ ■ ■ , Cfc)) 
t=0 

+ ^2(-^y +J ^(.Wi,(Tj],(T , ...,ai,...,Oj,.. .,a k ), 

i<j 

for uj G Sec(A k E*) and o"o, • • • , Cfc G Sec(r). d is a cohomology operator, that is, d 2 = 0. In 
particular, if /: M — > R is a real smooth function then df(a) = p(cr)f, for a € Sec(r). Locally, 

dx* = ^ Q e Q and de 7 = --C^e" A e^, 

where {e°} is the dual basis of {e a }. The above mentioned structure equations are but the 
relations d 2 x l = and d 2 e a = 0. We may also define the Lie derivative with respect to 
a section a of E as the operator d a : Sec(A k E*) — > Sec(A k E*) given by d a = i a o d + d o 
Along this paper, except otherwise stated, the symbol d stands for the exterior differential on a 
Lie algebroid. 



Morphisms 

Given a second Lie algebroid r': -E 1 ' — > M', a vector bundle map <!>: E ^ E' over M — > M' is 
said to be admissible if it maps admissible curves in E into admissible curves in E' , or equivalently 
if p' o $ = T99 o p. The map $ is said to be a morphism of Lie algebroids if Q*d9 = dQ*9 for 
every p-form 9 € Sec(A p £'*). Every morphism is an admissible map. 

In coordinates, a vector bundle map $>(x,y) = (ip 1 (x),^'g(x)y 13 ) is admissible if and only if 



k 



Moreover, such a map is a morphism if in addition to the above equation it satisfies 
p l — ^ - d — - - C ,c k <I> /3 $ 7 = 



Prolongation 

In what respect to Mechanics, the tangent bundle to a Lie algebroid, to its dual or to a more 
general fibration does not have an appropriate Lie algebroid structure. Instead one should use 
the so called prolongation bundle which has in every case the appropriate geometrical struc- 
tures [ZOOS]. 

Let (E, [ , ], p) be a Lie algebroid over a manifold M and v : P — > M be a fibration. For 
every point p £ P we consider the vector space 

T p E P = { (b, v) € E x x T P P I p(b) = T p v{v) } , 

where Tv : TP — » TM is the tangent map to z/ and ^(p) = x. The set T E P = \J p <zpT p E P 
has a natural vector bundle structure over P, the vector bundle projection t e being just the 
projection T E (b,v) = Tp(v). We will frequently use the redundant notation (p,b, v) to denote 
the element (6, v) € T p E P. In this way, the projection t e is just the projection onto the first 
factor. 
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The vector bundle r E : T E P -> P can be endowed with a Lie algebroid structure. The 
anchor map is the projection onto the third factor and will also be denoted by p, that is, the 
map p: 

T Ep Tp 

given by p(p, b, v) = v. To define the bracket on sections of T E P we will 
consider some special sections. A section Z £ Sec(T E P) is said to be projectable if there exists 
a section a £ Sec(E') such that Z(p) = (p, cr(V(p)), U(p)), for all p € P. Now, the bracket of two 
projectable sections Z\, Z2 given by Zi(p) = (p, (Ji{v(p)), Ui(p)), i = 1,2, is given by 

[Zi, Z 2 ](p) = Cp, K, a 2 ](i/(p)), [t/i, c/ 2 ](p)), p £ P. 

Since any section of T E P can be locally written as a C°°(M)-hnear combination of projectable 
sections, the definition of the Lie bracket for arbitrary sections of T E P follows. 

The Lie algebroid T E P is called the prolongation of v: P — > M with respect to E or the 
-E-tangent bundle to v. 

Given local coordinates (x l ,u A ) on P and a local basis {e a } of sections of E, we can define 
a local basis {X a , Va] of sections of T E P by 



X a (p) = [p,e a (u(p)),p l a ^J and V A (p) = (p»°'^a| )• 

If z = (p, b, v) is an element of T E P, with b = z a e a , then v is of the form v = p t OL z a -J~ + v A -^x, 
and we can write 

z = z a X a (p)+v A V A {p)- 

Vertical elements are linear combinations of {Va}- 

The anchor map p applied to a section Z of T E P with local expression Z = Z a X a + V A V A 
is the vector field on P whose coordinate expression is 

p(z) = p^4- + v a 







dx % ' du A 

The Lie brackets of the elements of the basis are given by 

[X a ,Xp]=Cr ap PC 1 , [X a ,V B ]=0 and [V A ,V B ]=0, 

and, therefore, the exterior differential is determined by 
dx i = PaX a , du A = V A , 

dX^ = --Cl R X a A X 13 , dV A = 0, 

where {X a , V A } is the dual basis to {<Yq,,V_a}. 



Prolongation of maps 

We consider now how to prolong maps between two fibrations v: P —* M and v' : P' — > M'. 
Let ^ : P — > P' be a map fibered over (p: M — > M'. We consider two Lie algebroids r: E — > M 
and t' : E' — > M' and a map $ : E — > E' fibered over ip. If $ is admissible, then we can define 
a vector bundle map T*^: T E P^T E 'P'hj means of 

T^( P ,b,v) = (*(p),*(&),rtt(«)). 

It follows that T*^ is also admissible. In [19] it was proved that T*^ is a morphism of Lie 
algebroids if and only if $ is a morphism of Lie algebroids. 

In particular, when E = E' and <I> = id we have that any map form P to P' fibered over the 
identity can be prolonged to a morphism T ld ^ which will be denoted simply by T^. We will 
also identify T E M (the prolongation of the 'fibration' id: M — > M with respect to E) with E 
itself by means of (m,b, p(b)) = b. With this convention, the projection onto the second factor 
of T E P is just T v: T E P -> E. It follows that Tv is a morphism of Lie algebroids. 
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3 Symplectic Mechanics on Lie algebroids 

By a symplectic structure on a vector bundle tt: F — ► M we mean a section to of A 2 F* which 
is regular at every point when it is considered as a bilinear form. By a symplectic structure on 
a Lie algebroid E we mean a symplectic section to of the vector bundle E which is moreover 
d-closed, that is duo = 0. A symplectic Lie algebroid is a pair (E,to) where E is a Lie algebroid 
and a; is a symplectic section on it. 

On a symplectic Lie algebroid (E, u>) we can define a dynamical system for every function 
on the base, as in the standard case of a tangent bundle. Given a function H £ C°°(M) there 
is a unique section au £ Sec(r) such that 

i aH (o = dH. 

The section an is said to be the Hamiltonian section defined by H and the vector field Xjj = 
p{(Jh) is said to be the Hamiltonian vector field defined by H. In this way we get the dynamical 
system x = Xh(x). 

A symplectic structure uona Lie algebroid E defines a Poisson bracket { , } on the base 
manifold M as follows. Given two functions F,G € C°°(M) we define the bracket 

{F,G} = uj(a F ,a G ). 

It is easy to see that the closure condition duo = implies that { , } is a Poisson structure on M. 
In other words, if we denote by A the inverse of to as bilinear form, then {F, G} = A(dF, dG). 
The Hamiltonian dynamical system associated to H can be written in terms of the Poisson 
bracket as x = {x,H}. 

Two important particular classes of symplectic dynamical systems on Lie algebroids are the 
following. 

Hamiltonian Mechanics (1 2 L 116] 

On T E E*, the ^-tangent to the dual bundle tt: E* -> M, we have a canonical symplectic 
structure. 

The Liouville section O € Sec((T E E*)*) is the 1-form given by 

(e,(fi,b, w )) = (fi,b). 

The canonical symplectic section Q, € Sec(A 2 (T E E*)*) is the differential of the Liouville section 

n = -de. 

Taking coordinates (x l ,/j, a ) on E* and denoting by {X a , V 13 } the associated local basis of sections 
T E E*, the Liouville and canonical symplectic sections are written as 

e = v a x a and n = x a a v a + \^cl p x a A X?, 

where {X a ,Vf3} is the dual basis of {X^V ^}. 

The Hamiltonian section defined by a function H G C°°(E*) are given in coordinates by 

and therefore, Hamilton equations are 

dx % j dH d[i a j dH 7 dH 

~dF =pa Wa ~dT = ~ Pa dx^~^ 

The Poisson bracket { , } defined by the canonical symplectic section $7 on E* is but the 
canonical Poisson bracket, which is known to exists on the dual of a Lie algebroid [5] and 
Hamilton equations thus coincide with those defined by Weinstein in [25] . 
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Lagrangian Mechanics 

The Lie algebroid approach to Lagrangian Mechanics builds on the geometrical structure of the 
prolongation of a Lie algebroid [H] (where one can develop a geometric symplectic treatment of 
Lagrangian systems parallel to J. Klein's formalism 

On the ^-tangent T E E to E itself we do not have a canonical symplectic structure. Instead, 
we have the following two canonical objects: the vertical endomorphism S: T E E — ► T E E which 
is defined by 

S(a,b,v) = (a,0,bl), 

where b^ a denotes the vertical lift to T a E of the element b G E, and the Liouville section, which 
is the vertical section corresponding to the Liouville vector field, 

A(o) = (a, OX). 

Given a Lagrangian function L G C°°(E) we define the Cartan 1-section 9 L G Sec((T E E)*) 
and the Cartan 2-section G Sec(A 2 (T E E)*) and the Lagrangian energy El G C°°(E) as 

6 L = S*{dL), Lu L = -d9 L and E L = £ A L - L. 

If the Cartan 2-section is regular, then it is a symplectic form on the Lie algebroid T E E, and we 
say that the Lagrangian L is regular. The Hamiltonian section Tl corresponding to the energy 
is the Euler-Lagrange section and the equations for the integral curves of the associated vector 
field are the Euler-Lagrange equations. 

If (x l ,y a ) are local fibered coordinates on E, (p l a ,C2p) are the corresponding local structure 
functions on E and {X a , V a } is the corresponding local basis of sections of T E E then SX a = V a 
and SV a = 0, and the Liouville section is A = y a V a . The energy has the expression El = 
-§^ky a — L, and the Cartan 2-section is 

d 2 L If d 2 L , d 2 L < 8L 



dy a dyP 



ya a yg , 1 ( 9 _lL J 9 L J + 9L \ vc 



AX? 



d 2 L 



from where we deduce that L is regular if and only if the matrix W a g = — — —, j is regular. In 

oy a oyP 

such case, the local expression of Tl is 

r L = y a x a + f a v a , 

where the functions f a satisfy the linear equations 

d y Pdy aJ ^ dx i dy a ^ y dy~f af3U Pa dx i 
Thus, the Euler-Lagrange equations for L are 

x ~ PaV ' dA dy- ) + <%7 C ^ y dx *-°- 

Finally, we mention that, as in the standard case, the relation between the Lagrangian and 
the Hamiltonian formalism is provided by the Legendre transformation Tl : E —> E* defined by 

(T L (a),b) = —L(a + tb)\ t=0 , 

for a,b G E with r(a) = r(b). Then it is easy to see that 

TT L *{Q) = 9 L and TT L *(tt) = oj l 

and therefore, in the regular case, the corresponding Hamiltonian sections are related by Tjj o 
Tl = TT l oTl. 
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4 Variational description 

While the Lie algebroid approach to geometric Mechanics builds on the geometrical structure of 
the T E E, it is well known that the origin of Lagrangian Mechanics is the calculus of variations. 
Integral curves of a standard Lagrangian system are those tangent lifts of curves on the base 
manifold which are extremal for the action functional defined on a space of paths. 

It is therefore interesting to find a variational description of Lagrange's equations for a Lag- 
rangian system defined on a more general Lie algebroid. The first steps in this direction where 
already done by A. Weinstein in [25] in the case of an integrable Lie algebroid (i.e. the Lie 
algebroid of a Lie groupoid) and by the author in |16|. [3], Finally, a formulation in the infinite 
dimensional manifold of curves was developed in [20j . 



The standard case 

Let us consider first the situation in the standard case of Lagrangian Mechanics, where E = TM. 
Given a Lagrangian function L 6 C°°{TM) we want to find those curves v : [to,t\] — > TM 
which are tangent prolongation of a curve in M, that is v = 7 for 7 = r o v, that connect 
to given points mo and mi in the base manifold M , and are extremal points of the action 
functional S(v) = J" 1 L(v(t))dt. One proceed as follows: given a solution v(t) = j(t) we 
consider variations v s (t) of v(t) such that vo(t) = v(t) and, for every fixed s, v s (t) = j s (t), 
where 7 s (i) = r(v s (t)). The infinitesimal variation is the vector Z(t) along v(t) given by 

Z(t) = 4zV s (t) , which obviously projects onto the vector field W(t) along j(t) given by 

s=0 

W(t) = 7357s (i) • The Euler-Lagrange equations 5L = are derived then by standard 

s=0 

manipulations of the condition for stationary points 



^S(v s ) = f\5L,W)dt. 
ds s=o J tQ 



Notice that 7 S determines v s , and hence W determines Z. This is clear in natural local 
coordinates (x^v*) on TM: we have that W = W^t)-^ and Z = W 1 -^ + W 1 -^. In classical 
notation 5x % = W l and 5v l = W l , which is but the well known rule for calculating the variation 
of the velocities as the derivative of the variation of the coordinates 

dv i = 6 (^-) =^(5x i ) = W i . 
\dt J dV 1 

Finally notice that due to the fixed endpoints condition we have that W(to) = W{t\) = 0. 

Geometrically, things are a bit more difficult. The vector field W(t) is a vector field along 
the curve 7(i), and hence, it is a curve in TM over j(t). If we take the tangent lift, the curve W 
is a curve in TTM over Wit) and therefore defines a vector field along W(t), instead of a vector 
field along v{t) = j(t). Therefore, the variation vector field Z(t) is not just W(t) since they 
are defined at different points. A further operation is needed and this is the so called canonical 
involution or Tulczyjew involution. It is a map xtm ■ TTM — > TTM such that 

/ d 2 (3 , \ d 2 (5 , . 

for every map 13: R 2 M locally defined in a neighborhood of the origin. It follows that the 
variation vector field Z is not W but it is 



Z = X tm(W). 
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In many situations this kind of variations is obtained in terms of the flows of vector fields. 
Given a vector field X on the manifold M we consider its flow {ips} and then we define a variation 
of v(t) by v s (t) = Tip s (v(t)). It is clear that they are admissible variations, being the base 
variations 7 s (i) = ipsfait)), that W(t) = X(7(i)) and Z(t) = X c (v(t)). In this expression, 
X c S 3C(TM) is the complete or tangent lift of X, whose flow is {Tip s } and which can be 
defined in terms of the canonical involution by means of 

X c (v) = X tm{TX{v)), for all v G TM. 

Using this kind of variations, the Euler-Lagrange equations can be easily found to be 

— (6 L ,X C )= C x cL, 
where £ denotes the Lie derivative. 



The general case 

In the general case of a Lagrangian system on an arbitrary Lie algebroid E we can follow a similar 
path. We consider a Lagrangian L G C°°(E) and the action 

5(a) = / 1 L{a(t))dt 

defined on the set of admissible curves on E with fixed base endpoints tuq at to an d mi at t±. 
We look for a variational principle for the Euler-Lagrange equations, that is we have to specify 
boundary conditions and a class of variations such that the critical points of the action are 
precisely those curves satisfying Lagrange equations. As we will see such variations are related 
to complete lifts of sections of E. 

Every section i] of E can be naturally lifted to a section of T E E in two different ways: the 
vertical lift if and the complete lift rf . The structure of Lie algebroid in T E E is determined 
by the brackets of such sections, 

[ V c , a c ] = [a, r/] c , [r/ c , a v ] = [a, ^ and [^ v , a v ] = 0. 

This relations were used in [15] to define the Lie algebroid structure, so that we mimic (and 
hence extend) the properties of complete and vertical lifts in the tangent bundle, which are on 
the base for the geometric formalism in the calculus of variations. 

In local coordinates, if rj = r/ a e a is a local section of E then the vector field associated to its 
complete lift has the local expression 

and the associated vector field has the expression 

where / = Pay a §§i- More generally, one can define the complete lift of a time-dependent section, 
which has a similar expression as long as one defines / = §{" + PaV a §^- 

Using the properties of complete and vertical lifts, it is easy to see that the Euler-Lagrange 
equations ir L ^L = dE^ can also be written in the form 



dr L {0 L ,cr c ) = d a cL, 
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for every time-dependent section a of E. From this expression one can deduce that the infinite- 
simal variations that one must consider are precisely the vector field associated to the complete 
lifts of sections of E. 

The above observation is not only a formal statement, but can be carried on precisely in 
terms of the flow associated to a time-dependent section (see [21 02]). For simplicity in the 
exposition I will consider only time-independent sections. If rj is a section of E then the flow $ s 
of the vector field p(rj c ) S 3L{E) projects to the flow ip s of the vector field p(rj) E X(Af). For 
every fixed s, the map <J> S is a vector bundle map which is a morphism of Lie algebroids over (p s . 
The pair (<J> S , <p s ) is said to be the flow of the section rj £ Sec(E), and we have that 

de = ^0 

as s=o 

for every tensor field 9 over E. 

Given an admissible curve a(t) we consider a section a of E and its flow (<& s , ip s ) the variations 
a s (t) = $ s (a(t)), which are also admissible curves, since $ s are morphisms of Lie algebroids. If 
we moreover consider sections rj vanishing at the endpoints, 77(7710) = 77(7711) = 0, then the varied 
curve has fixed endpoints, r(a s (to)) = niQ and r(a s (ti)) = niQ. 

Notice that, in general, there are more general variations preserving the admissibility of curves 
than those considered here. Nevertheless, we have to chose exactly the ones we have chosen: if 
we consider a restricted class of variations, we will get unspecified dynamics and if we consider 
a more general class of variations we will get some constraints. One can clearly see this fact in 
the case of a Lie algebra, where every curve is admissible and hence every variation preserves 
admissible curves. 



The canonical involution 

In the argument given above, in order to define a variation we need a section of E defined 
in a neighborhood of the base path. As in the case of the standard Lagrangian mechanics, 
a different procedure consists in using the canonical involution for defining variations. 

Indeed, the canonical involution can also be defined on any Lie algebroid E (see |12] for the 
details). That is, there exists a canonical map xe'- T e E — » T E E such that x\ = ^ and it is 
defined by X.e(a, b, v) = (b,a,v), for every (a,b,v) € T E E, where v € T^E is the vector which 
projects to p(a) and satisfies 

v§ = v9 + d6(a,b) 

for every section 9 of E*, where 9 G C°°(E) is the linear function associated to 9. In local 
coordinates the canonical involution is given by 

XE(x\y a ,z a ,v a ) = (x\z a ,y a ,v a + C^z /3 y"'). 

From this expression is clear that the complete lift of a section 7] £ Sec(E) can be given in terms 
of the canonical involution by 

rf(a) = XE\Tr](aj\ for all a G E. 

This formula suggests to consider the following map. Given an admissible curve a: R — > E 
over 7 = r o a we consider the map S a from sections of E along 7 to sections of TE along a, i.e. 
E a : Sec 7 (£) -> Sec a (T£'), given by 

^a(o-) = p 1 (xE(cr,a,a)). 

The local expression of the map H a is 







a(a)(t) = pMt))a a (t)^- m + (a a (t) + C^(t))a^ ' {t)a> (t) 



a(t) V ^ " " >)dy 







a(t) 



where a and a have the local expression a(t) = (7* (t) , a a (t)) and a(t) = (7*(£), a a (t)). 



10 



E. Martinez 



5 The manifold of .E-paths 

To get some more insight into the variational principle that we have obtained, we can analyze 
the situation from the point of view of the geometry of the infinite dimensional manifold of 
admissible curves. 

Homotopy of .E-paths 

Let I = [0, 1] and J = [io,ii], and denote the coordinates in M 2 by (s,t). Given a vector bundle 
map TM 2 — » E, denote a(s,t) = $>(dt\( s ,t)) an d b(s,t) = $(9 s |( S)t O, so that we can write 
$ = adt + bds. 

Definition 1. Two E-paths ao and a\ are said to be E-homotopic if there exists a morphism 
of Lie algebroids $ : TI x TJ — > E, ® = adt + bds, such that 

a(0,t) =a (t), b(s,to)=0, 
a(l,t) = ax{t), b(s,h) = 0. 

We will say that $ is an E-homotopy from the E-path ao to the E-path a\. 

It follows that the base map is a homotopy (in the usual sense) with fixed endpoints between 
the base paths. Notice that a(s,t) is a variation of a(0, t) and one should think of b(s,t) as the 
vector generating the variation. 

Theorem 1 ([ID]). The set of E-paths 




is a Banach submanifold of the Banach manifold of C 1 -paths whose base path is C 2 . Every 
E-homotopy class is a smooth Banach manifold and the partition into equivalence classes is 
a smooth foliation. The distribution tangent to that foliation is given by a £ A(J, E) i— > F a 
where 



and the codimension of F is equal to dim(E) . The E-homotopy equivalence relation is regular 
if and only if the Lie algebroid is integrable (i.e. it is the Lie algebroid of a Lie groupoid). 

The space of E-paths 

Therefore, on the same set A( J, E) there are two natural differential manifold structures: as 
a submanifold of the set of C 1 paths in E, which will be denoted just A(J, E), and the structure 
induced by the foliation into E-homotopy classes, which will be denoted V{J,E). We will refer 
to it as the space of E-paths on the Lie algebroid E. The structure of A( J, E) is relevant when 
one wants to study the relation between neighbor E-homotopy classes, as it is the case in the 
problem of integrability of Lie algebroids to Lie groupoids. The structure of V(J, E) is just the 
structure that one needs in Mechanics, where one does not have the possibility to jump from 
one E-homotopy class to another. 

Notice that every homotopy class is a connected component of V(J,E), and the identity 
defines a smooth map i: V(J,E) — > A(J,E) which is an (invertible) injective immersion. The 
image by i of a leaf is an immersed (in general not embedded) submanifold of A(J,E). The 
tangent space to V(J,E) at a is T a V{J,E) = F a . The topology of V{J,E) is finer than the 
topology of A(J, E). In particular, if G : A(J, E) — > Y is a smooth map, then Goi: V(J, E) Y 
is also smooth. 





F a = { E a (a) G T a A(J, E) \ a{t ) = and a(h) = } . 
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Variational description 

With the manifold structure that we have previously defined on the space of £7-paths, we can 
formulate the variational principle in a standard way. Let us fix two points mo, mi 6 M and 
consider the set V(J,E)™* of those E'-paths with fixed base endpoints equal to mo and m±, 
that is 

V(J,E)™1 ={aEV(J,E) | r(a(to)) = mo and r(a(ii)) = m ± } . 

We remark that V(J, E)™* is a Banach submanifold of V(J,E), since it is a disjoint union of 
Banach submanifolds (the .E-homotopy classes of curves with base path connecting such points). 
On the contrary, there is no guaranty that the analog set A(J, E)™* is a manifold (see |22|). 

Theorem 2 (|20j). Let L £ C°°(E) be a Lagrangian function on the Lie algebroid E and fix 
two points m,Q,mi £ M. Consider the action functional S: V(J,E) — > ]R given by S(a) = 
L{a{t))dt. The critical points of S on the Banach manifold V{ J, are precisely those 

elements of that space which satisfy Lagrange 's equations. 

Lagrange Multipliers 

We can also analyze the problem by using Lagrange multipliers method by imposing a condition 
on A( J, E) which represents the constraint that our S-paths are in a given .E-homotopy class. 
This is connected with the theory of Lin constraints [1] . 

We consider only the case of an integrable Lie algebroid, since in the contrary we will not 
have a differential manifold structure in the set of E-homotopy equivalence classes. In this 
case, the foliation defined by the E'-homotopy equivalence relation is a regular foliation so that 
quotient G = A(J,E)/ ~ has the structure of quotient manifold and the quotient projection 
q: A(J,E) — » G is a submersion. Defining the source and target maps by s([a]) = r(a(io)) and 
i([a]) = r(a(ii)), the unit map e: M — > G by e(m) = [0 m ], where m denotes the constant curve 
with value £ E m , and the multiplication induced by concatenation of S-paths, we have that 
G is the source simply-connected Lie groupoid with Lie algebroid E. See [10] for the details. 

Given g 6 G, we can select the curves in an .E'-homotopy class as the set q _1 {g)- Therefore we 
look for the critical points of the functional S{a) = L(a(t)) dt defined in A{J, E), constrained 
by the condition q(a) = g. Since q is a submersion, there are not singular curves for the constraint 
map, and we can use Lagrange multiplier method in the standard form pQ. 

Theorem 3 ([20]). Let S : A(J,E) R, be the action functional S(a) = J"* 1 L(a(t)) dt. An 
admissible curve a € A(J,E) is a solution of Lagrange's equations if and only if there exists 
fi £ T*G such that dS(a) = /i o T a q. 

A more convenient setting for this constrained problem consists in fixing one of the endpoints. 
Given mo € M, the subset A{ J, E) mQ of those E-paths whose base path start at mo, A{ J, E) mo = 
{ a G A( J, E) | r(a(to)) = mo }, is a smooth Banach submanifold of A(J, E). On A( J, E) mo we 
define the map p: A(J,E) mo — > s~ 1 (mi) by p(a) = L g -i(q(a)). With the help of this map, the 
constraint reads p(a) = e{m\), because an £"-path is in q~ 1 (g) if and only if it is in p~ l (e{m{)). 
Then the tangent space to A( J, E) mo at a S A(J, E) mo is 

T a A(J, E) mo = { E a (a) | a{t ) = } . 

p is a submersion and the tangent map T a p: T a A(J, E) mo — > E mi to p at a point a S p _1 (e(mi)), 
is given by the endpoint mapping 

T a p(E a (cr)) = a(tt) 

for every a € Sec^(E) such that o~(to) = 0. If we now apply Lagrange multiplier theorem we 
obtain the following result. 
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Theorem 4 ([20]). Let S mo be the restriction of the action functional to the submanifold 
A(J,E) mo . An admissible curve a £ A(J,E) mo is a solution of Lagrange's equations if and 
only if it there exists A 6 E^ such that dS mo (a) = A o T a q. The multiplier A is given explicitly 
by\ = e L (a(h)). 

6 Morphisms and reduction 

One important advantage of dealing with Lagrangian systems evolving on Lie algebroids is that 
the reduction procedure can be naturally handled by considering morphisms of Lie algebroids, 
as it was already observed by Weinstein [25] . We study in this section the transformation laws 
of the different geometric objects in our theory and we apply these results to the study of the 
reduction theory. 

Mappings induced by morphisms 

We recall that admissible maps are precisely those maps which transforms admissible curves into 
admissible curves. Therefore an admissible map <£: E — > E' induces a map between E-paihs by 
composition a \—> $oa. We prove now that such a map is smooth provided that $ is a morphism. 

More precisely, let <3>: E — > £' be an admissible map. It is easy to see that, $ is a Lie 
algebroid morphism if and only if T<£> o H a (<r) = S$ a(^ ° c) for every E-p&th a and every 
section a along the base path r o a. It follows that morphisms transform vectors tangent to the 
foliation into vectors tangent to the foliation, and hence they induce a smooth map between 
path spaces. 

Proposition 1 ([20]). Given a morphism of Lie algebroids <I>: E — > E' the induced map 
<3? : V(J,E) — > V(J, E') given by &(a) = <I> o a is smooth. Moreover, 

• If & is fiberwise surjective then $ is a submersion. 

• If & is fiberwise injective then & is a immersion. 

As a consequence, the variational structure of the problem is not broken by reduction. On 
the contrary, reduction being a morphism of Lie algebroids, preserves such structure. The 
above results says that morphisms transforms admissible variations into admissible variations. 
Therefore, a morphism induces relations between critical points of functions defined on path 
spaces, in particular between the solution of Lagrange's equations. 

Reduction of the variational principle 

Consider a morphism $ : E — > E' of Lie algebroids and the induced map between the spaces of 
paths : V{ J, E) — > V(J, E'). Consider a Lagrangian L on E and a Lagrangian V on E' which 
are related by <&, that is, L = L' o&. Then the associated action functionals S on V(J, E) and S' 
onV(J,E') are related by <l, that is S' o l> = S. Indeed, 

5'($(a)) = 5'($oa)= [ o $ o a)(t) dt = [ ' (L o a){t) dt = 5(a). 

Jt Jto 

The following result is already in [25] but the proof is different. 

Theorem 5 (\25\ I20j). Let $: E — > E' be a morphism of Lie algebroids. Consider a Lagran- 
gian L on E and a Lagrangian L' on E' such that L = L' o $. If a is an E-path and a' = $ o a 
is a solution of Lagrange 's equations for L' then a itself is a solution of Lagrange 's equations 
for L. 
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Proof. Since S' o <l = S we have that ( dS'($(a)) , T a $(v) ) = (dS(a),v) for every d € 
T a V(J,E)™*. If $(a) is a solution of Lagrange's equations for L' then dS'($>(a)) = 0, from 
where it follows that dS(a) =0. ■ 

From the above relations between the action functionals it readily follows a reduction theorem. 

Theorem 6 (|20j). Let <3>: E — > E' be a fiberwise surjective morphism of Lie algebroids. Con- 
sider a Lagrangian L on E and a Lagrangian L 1 on E' such that L = L' o If a is a solution 
of Lagrange 's equations for L then a 1 = $ o a is a solution of Lagrange 's equations for L' . 

Proof. Since S' o $ = S we have that (dS'($(a)) ,T *(«) } = (dS(a),v) for every u 6 
T a V(J,E)™i. If $ is fiberwise surjective, then $ is a submersion, from where it follows that <3? 
maps critical points of S into critical points of S', i.e. solutions of Lagrange's equations for L 
into solutions of Lagrange's equations for L' . ■ 

Reduction of the symplectic form and the dynamics 

Reduction can also be studied in the context of the symplectic formalism on Lie algebroids, 
see [8] or [9] for the details. 

Proposition 2. Let $ : E — > E' be a morphism of Lie algebroids, and consider the <&-tangent 
prolongation o/ i.e 7~*<I>: T E E — > T E E' . Let S and S' , and A and A', be the the vertical 
endomorphisms and the Liouville sections on E and E' , respectively. Then, 

o A = A' o $ and oS = S'o T®&. 

Proposition 3. Let L € C°°(E) be a Lagrangian function, 8l the Cartan form and u>l = —dO^. 
Let E — > E' be a Lie algebroid morphism and suppose that L = L' o with L' 6 C°°(E') 
a Lagrangian function. Then, we have 

(T*$)*6 L , = 9 L , {T^^Yuj v =uj l and (T*$)*E L , = E L . 

The transformation of the symplectic equation is easily found by means of standard argu- 
ments, and we find 

{T*<S>y{i TL ,u v - dE v ) - (i TL uj L - dE L ) = uj l ,{T l , o $ - o T L ,T*$( • )), 

It follows that, if $ is a fiberwise surjective morphism and L is a regular Lagrangian on E, 
then L' is a regular Lagrangian on E' (note that : T E E -» T E ' E' is a fiberwise surjective 
morphism) we have that the dynamics of both systems is uniquely defined, and it is related as 
follows. 

Theorem 7. Let Tl and be the solutions of the dynamics defined by the Lagrangians L 
and L' , respectively, with L = L' o $ . If $ is a fiberwise surjective morphism and L is a regular 
Lagrangian, then L' is also a regular Lagrangian and 

r $ $or L = T v 0$. 

Finally, by introducing some constraints one can also study nonholonomic mechanical systems 
on Lie algebroids. See [U [9] for the general theory and results on reduction for nonholonomic 
systems. 
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Examples 

We present here some examples where the reduction process indicated above can be applied. 
See [2U1 [SJ E] for more examples. 

Lie groups. Consider a Lie group G and its Lie algebra q. The map <E> : TG — > g given by 
&(g,g) = g~ l g is a morphism of Lie algebroids, which is fiberwise bijective. As a consequence 
if L is a left-invariant Lagrangian function on TG and L' is the projected Lagrangian on the Lie 
algebra q, that is L(g,g) = L'(g~ l g), then every solution of Lagrange's equations for L projects 
by $ to a solution of Lagrange's equations for V . Moreover, since $ is surjective every solution 
can be found in this way: if the projection = g(i) -1 <?(i) of an admissible curve (g(t),g(t)) 
is a solution of L', then (g(t),g(t)) is a solution for L. Thus, the Euler-Lagrange equations on 
the group reduce to the Euler-Poincare equations on the Lie algebra, both being symplectic 
equations in the Lie algebroid sense. 

Lie groupoids. Consider a Lie groupoid G over M with source s and target i, and with 
Lie algebroid E. Denote by T S G — > G the kernel of Ts with the structure of Lie algebroid 
as integrable subbundle of TG. Then the map $ : T S G — > E given by left translation to 
the identity, <&(v g ) = TL g -i(v g ) is a morphism of Lie algebroids, which is moreover fiberwise 
surjective. As a consequence, if L is a Lagrangian function on E and L is the associated left 
invariant Lagrangian on T S G, then the solutions of Lagrange's equations for L project by <E> 
to solutions of the Lagrange's equations. Since $ is moreover surjective, every solution can be 
found in this way. 

Group actions. We consider a Lie group G acting free and properly on a manifold Q, so 
that the quotient map tt : Q — > M is a principal bundle. We consider the standard Lie algebroid 
structure on E = TQ and the associated Atiyah algebroid E' = TQ/G — > M. The quotient 
map $: E — > E', &(v) = [v] is a Lie algebroid morphism and it is fiberwise bijective. Every 
G- invariant Lagrangian on TQ defines uniquely a Lagrangian V on E' such that V o $ = L. 
Therefore every solution of the G- invariant Lagrangian on TQ projects to a solution of the 
reduced Lagrangian on TQ/G, and every solution on the reduced space can be obtained in 
this way. Thus, the Euler-Lagrange equations on the principal bundle reduce to the Lagrange- 
Poincare equations on the Atiyah algebroid, both being symplectic equations in the Lie algebroid 
sense. 

7 Optimal control theory 

As it is well known, optimal control theory is a generalization of classical mechanics. It is 
therefore natural to see whether our results can be extended to this more general context. The 
central result in the theory of optimal control systems is Pontryagin maximum principle. The 
reduction of optimal control problems can be performed within the framework of Lie algebroids, 
see [17] , This was done as in the case of classical mechanics, by introducing a general principle 
for any Lie algebroid and later studying the behavior under morphisms of Lie algebroids. 

Pontryagin maximum principle |17] 

By a control system on a Lie algebroid r: E — > M with control space tt: B — > M we mean 
a section a of E along tt. A trajectory of the system a is an integral curve of the vector 
field p(o) along tt. Given an index function L G C°°(B) we want to minimize the integral of L 
over some set of trajectories of the system which satisfies some boundary conditions. Then 
we define the Hamiltonian function H £ C°°(E* Xm B) by H(fj,,u) = (n,a{u)) — L{u) and 
the associated Hamiltonian control system an (a section of T E E* along pr x : E* Xm B — ► E*) 
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defined on a subset of the manifold E* x m B, by means of the symplectic equation 
r tt = dH. 



(*) 



The integral curves of the vector field p((Jh) are said to be the critical trajectories. 

In the above expression, the meaning of i fTH is as follows. Let <I> : E — > E' be a morphism 
over a map ip: M — » M' and let 77 be a section of E' along p. If u is a section of /\ p E'* then 
Zjju; is the section of /\ p l E* given by 

(i v u) m (ai, Op-x) = uJ v(n) (r](m), $(a x ), . . . , $(a p _i)) 

for every m G M and ai, . . . , a p _i G E m . In our case, the map $ is Tpr 1 : T E (E* Xjf 1?) — > 
T E E* , the prolongation of the map pr x : E* Xm B —> E* (this last map fibered over the identity 
in M), and <7# is a section along pr x . Therefore, v H fJ — (Iff is a section of the dual bundle to 
T E (E* x M B). 

It is easy to see that the symplectic equation (j*j) has a unique solution defined on the following 
subset 

S H = {(fi,u) G E* x M B \ (dH(fi,u),V) = for all V G kerTp^}. 

Therefore, it is necessary to perform a stabilization constraint algorithm to find the integral 
curves of an which are tangent to the constraint submanifold. 

In local coordinates, the solution to the above symplectic equation is 



oh 



8H 

dp a 



dH 

'dx 1 



7 OH 



defined on the subset where 
dH 

0, 



du A 



and therefore the critical trajectories are the solution of the differential-algebraic equations 
.• dH 



l dpa 



: 



dH 

'dx 1 



dH 

dfi/3 



dH 

du A 



Notice that M- = a a . 

One can easily see that whenever it is possible to write [i a = Pip l a then the above differential 
equations reduce to the critical equations for the control system Y = p(a) on TM and the 
index L. Nevertheless it is not warranted that p is of that form. For instance in the case of a Lie 
algebra, the anchor vanishes, p = 0, so that the factorization p a = Pip l a will not be possible in 
general. 



Reduction 

Consider two optimal control systems, with data (B,E,M,a,L) and {B' , E' , M' ,a' , L'). Let 
$ : E — > E' be a fiberwise bijective morphism of Lie algebroids over a map <p. Then <£ m : E m — > 
E r f m \ is invertible for every m G M, and we can consider the contragredient map <3? c : E* — > 

E*. This is the vector bundle map over p whose restriction to the fiber over m G M is 
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$m = ®m 1 -E* m ^ E'* {m) , given by ( $ c (/x) , a ) = (n,$J(a)) for every p € E* m and every 

ip(m) 

If we further have a fibered map tp: B — > B' over the same map (p, then we can define the 
map ^: E* x M B -> E'* x M > B' by 

«) = ($ c (/i), (/u, «) eE*x M B, 

and we have the following transformation properties. 

Proposition 4. We have the following properties 

1. (T*$ c )*9' = e. 
2 (T*$ c )*0' = ^. 
3. If V o if) = L, then H'o^ = H. 

Proof. For every (H,b,v) € T E E* we have 

( (r** c )*e' , (m, b, u) > = ( e' , r*$ c (/i, b, «) ) = ( e' , $(b), r$ c w) ) 

= (<I> c (/,),<i>(b)} = (^,b} = (e,(/,,b,t;)}, 

which proves the first. 

The proof of the third is similar and will be omitted. For the second just take into account 
that 7~ $< I> C is a morphism. ■ 

By mean of standard manipulations one can easily show that 

from where the following theorem, which establishes the relation between critical trajectories of 
the two related optimal control problems, readily follows. 

Theorem 8. Let ip: B — > B' and <E>: E — > E' be fibered maps over the same map cp: M — > M' , 
and assume that tjj is fiberwise submersive and $ is a morphism of Lie algebroids which is 
fiberwise bijective. Let L be an index function on B' and L' be an index function on B' such 
that L = L' o ip and let an and a^ 1 the corresponding critical sections. Then we have that 
^(S H ) C S h > and 

T*$ c o a H = oh' ° * 

on the subset Sh- 

As a consequence, the image under ^ of any critical trajectory for the index L is a critical 
trajectory for the index V . 

As an application of the above result we can consider the case of reduction by a symmetry 
group (with a free and proper action). Indeed, applying this result to E = TQ, B = B, M = Q 
and E' = TQ/G, B' = B/G, M' = Q/G, with ^(b) = [6], $(y) = [v], <p(q) = [q], the quotient 
maps and index L, Z/([b]) = L(b) (so that L = L' o ip) we have that any critical trajectory for L 
in Q is a critical trajectory for L' in the reduced space M'. 

Finally notice that the case of Hamiltonian mechanics corresponds to B = E and a = id, and 
hence the set Sh is the graph of the Legendre transform. Therefore, the above results about 
optimal control produce also results about reduction of Hamiltonian systems. It would be nice 
to have a similar result for fiberwise surjective (no necessarily fiberwise invertible) morphisms. 
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